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I. INTRODUCTION 



At high temperatures the electroweak symmetry is restored. Since the baryon violating processes are unsuppressed 
at high temperatures, the observed baryon asymmetry of the universe has finally been determined at the electroweak 
phase transition 

In recent years quantitative studies of the electroweak phase transition have been carried out by means of resummed 
perturbation theory and lattice Monte Carlo simulations |2j- In the SU(2)-Higgs model for Higgs masses (to//) 
below 50 GeV, the phase transition is predicted by the perturbation theory to be of first order. However, it is difficult 
to give a definite perturbative statement for physically more interesting masses, e.g. mu > 80 GeV. Due to the bad 
infrared properties of the theory, the perturbative approach breaks down in this parameter region. A systematic and 
fully controllable treatment is necessary, which can be achieved by lattice simulations. 

For smaller Higgs boson masses (to// < 50 GeV) the phase transition is quite strong and relatively easy to study on 
the lattice. For larger to// (e.g. to// = 80 GeV) the phase transition gets weaker, the lowest excitations have masses 
small compared to the temperature, T. From this feature one expects that a finite temperature simulation on an 
isotropic lattice would need several hundred lattice points in the spatial directions even for Lt = 2 temporal extension. 
These kinds of lattice sizes are out of the scope of the present numerical resources. 

One possibility to solve the problem of these different scales is to integrate out the heavy, 0{T) modes perturbatively, 
and analyse the obtained theory on the lattice. This strategy turned out to be quite successful, and both its perturbative 
and lattice features have been studied by several groups j^- ||ll[]. Even more: predicts that somewhere above 80 
GeV Higgs mass the first order phase transition does not take place any further, the two phases can be continuously 
connected. Ref. [l9|] gives estimates of the end-point of the phase transition in the framework of the reduced 3-d 
approach. 

With this paper we follow another approach (analytic and Monte Carlo) to handle this two-scale problem. We will 
use the simple idea that finite temperature field theory can be conveniently studied on asymmetric lattices, i.e. lattices 
with different spacings in temporal (at) and spatial (a^) directions. This method solves the two-scale problem in a 
natural way ||20|] . Another advantage is, well-known and used in QCD, that this formulation makes an independent 
variation of the temperature (T) and volume (V) possible. The perturbative corrections to the coupling anisotropics 
are known in QCD (see refs. |2|,|2^). Performing a similar analysis for the SU(2)-Higgs model, we have presented 
in our earlier paper | p3[ the perturbative corrections to the coupling anisotropics for this case, too. Here we give 
details of the perturbative calculation of |p3|| and by numerical simulation on lattices with anisotropic lattice spacings 
calculate the coupling asymmetries for the practically reasonable parameters of the SU(2)-Higgs model. As we will 
show at TO// ~ 80 GeV the one-loop perturbative and the non-perturbative coupling asymmetries agree very well. 

There is an essential difference between pure gauge theories and the SU(2)-Higgs model. In the former case any 
value (within a certain range) of the space and time coupling constant ratio does correspond to a meaningful theory, 
the actual value of the ratio corresponds to a definite value of the ratio of the space and time lattice spacings. On the 
other hand, in case of the SU(2)-Higgs model for a fixed value of the space and time gauge coupling ratio (determining 
the ratio of the space and time lattice spacings) one has to fix the ratio of the space and time hopping parameters to 
a definite value, in order to ensure that the theory makes sense. A convenient way to do this is to require that the 
ratio of space and time gauge boson masses should be equal to the ratio of space and time Higgs boson masses. Such 
a choice of the parameters is a precondition to both the perturbative calculations and the numerical simulations. 

The plan of this paper is as follows. Section II. deals with the perturbative analysis. In subsection II. A we give 
the lattice action of the model on asymmetric lattices and discuss perturbation theory in the anisotropic lattice case. 
Subsection II. B contains the calculation of the critical hopping parameter and of the wave function quantum correction 
terms, which give the quantum corrections to the anisotropy parameters. In subsection II. C a discussion of the finite 
temperature continuum limit is given. The optimal choice of the ratio of space and time lattice spacings is determined 
by perturbative techniques. Section III. contains our non-perturbative analysis. Subsection HI. A gives the basic 
points of our MC simulations. Subsection III.B deals with the mass determinations from the correlation functions. In 
subsection III.C we present the results on Wilson loop simulations and the static potential. In subsection III.D we 
finally evaluate the non-perturbative asymmetries and compare them with the perturbative results. Section IV. is a 
summary and outlook. 

II. LATTICE ACTION AND PERTURBATION THEORY 

In this section we discuss lattice perturbation theory. Since we did not find the Feynman rules for the anisotropic 
lattice spacing case in the literature, we present some details of lattice perturbation theory. After that we determine 
the critical hopping parameter and the anisotropics in one-loop perturbation theory. The continuum limit and the 
optimal choice of the ratio of space- and time-like lattice spacings is also discussed. 
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A. Action in continuum notation, gauge fixing, propagators 



For simplicity, we use equal lattice spacings in the three spatial directions (a.; — as, i — 1,2, 3) and another spacing 
in the temporal direction (04 — at). The asymmetry of the lattice spacings is characterized by the asymmetry factor 
^ = as I at- The different lattice spacings can be ensured by different coupling strengths in the action for time- like and 
space-like directions. The action reads 



S[U, = /3, ^ M - ^Tr Usp\ + /5t E (1 - ^^r Utp 

Sp ^ ^ tp ^ 

+ EUTr(v.+ ^,) + A" 



-Ks Tr {(p^^f^U.^,^, ifx) - KtTr {'^^^-^U^^^i (p^) \ , (1) 
^1=1 J 

where A stands for the lattice points, JT^^^ denotes the SU(2) gauge link variable, Usp and Utp the path-ordered product 
of the four Ux,i_i around a space-space and space-time plaquette, respectively. The symbol ipx stands for the Higgs 
field. 

The values of anisotropics defined as 

2 A 2 '^t 

are choosen to correspond to given values of the asymmetry ^. In perturbation theory this can be ensured order by 
order in the loop expansion, requiring that in the limit as, at — » with the ratio ^ = as /at fixed, certain physical 
quantities show rotation symmetry on submanifolds of the bare coupling space satisfying 7^ — const, 7^ — const. 
This procedure leads to a formal double expansion in and A of the anisotropics: 

7^ = [1 + cpiOg^ + b^iOX + 0{a\ A^)] , il = e[l + + 6«(0A + 0{a\ A^)] . (3) 

Here g is the bare gauge coupling of the theory with symmetric lattice spacings in standard notation. (Note that 
C(}{1) = Ck(1) — 6/3(1) — 6^(1) — 0.) In this double expansion we use the formal power counting A g^. In general, 
fixing 7/3 (^) and 7k(^) to ensure rotation symmetry should be done non-perturbatively. In the non-perturbative 
framework the definition of ^ is given as the ratio of space and time direction lattice unit correlation lengths. This 
non-perturbative analysis will be the topic of section III., where we choose the values of the bare copuling ratios (^ 
to ensure that the Higgs and gauge boson correlation lengths in physical units be the same in the different directions. 
This idea can be applied in perturbation theory as well (see e.g. [p2|), and we will follow this method in our analysis, 
too. 

Elaborating perturbation theory we follow the usual steps (see e.g. |24|, ||2^ for the isotropic SU(2)-Higgs model). 
The only complication is that we have to keep track of the different lattice spacings and couplings. 
First we consider the gauge part of the action. We will use the same notation as applied by the calculation for 
the pure gauge theory, 

Ux,f, = exp (ia^.gJ^A'^ix)^ , (4) 

where r is summed over 1,2,3, while /i is not summed, moreover a^ — as, g^, — gs for /i = 1, 2, 3 and — at, gi — gt. 
We have also 

C5s(0 5t(0 
as the connection to the lattice parameters. The expansions for gs, gt read: 

9s{£.f - g\i csiOg' bsiOX + 0{g\ X'j), (6) 



gtiO' - .9'(i + ctiOg' + hm + 0{g\ X')), (7) 



where g^ = gU^ = 1) = = 1) and cp{0 = ct{0 - c.(C), bpiO = hiO - bs{0- 
We write 



Ux,^^al{x)+iTral{x), (8) 



3 



where 



with = ^ A^^{x) A'^^{x) . The expansion is given by 



= 1 - ^-^^A:AI + ^-^^AIAIAIK + 0{g^), (10) 



alix) = ^Al - I^^A^A^A^ + 0{g% (11) 
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Inserting eq. (|^) into the plaquette parts of the lattice action we get the parts of the action containing odd numbers 
of gauge boson fields [S°f'^) and even numbers of gauge boson fields (5'p['™). They read: 



Spi"^ = - X! X! ^^^'^ {eprs(a°(a; + fia^,)aF^{x) + al^{x)al{x + iia^))aJ'^{x + va^)al{x) 

-eprsaP^j,{x)al{x + /ia^)(a°(x)aj;(a; + va^) + a°(x + i'ay)al{x))] , (12) 

geven ^ _ ^ ^ {a°^{x)al{x + (j.a^,)a°^{x + i>a^)a°(a;) - a°(a;)a"(a; + Aia^)a^(a; + Oa^,)al{x) 

-a°(x + i>a^)a°(a;)aj^(a;)a'j^(a; + fia^) + aj^(a;)a'j^(a; + fiap)al{x + i>a^)a^(a::) 
+ (a°(a; + jla^_,)a''^{x) + a°^{x)al{x + fla^_,)){al(x)a''^{x + va^) + a°{x + C'ai,)al{x)) 
+ i-6ps6rt + SptSrs)a^{x)al{x + fia^)a''^{x + va^)a!'^{x)] , 

(13) 

where p,r,s,t=l,. . . , 3 and Pfj_i, is equal to for space indices and equal to Pt for one space and one time index. 
The integration measure for the gauge variables also contributes to the action. 

d'U,,^. = ^d'a'^ix)d{al{x)alix) - 1) 
where capital letters run from to 3 and a sum over T=0,. . . , 3 is understood. The contribution to the action reads: 



^™ — E E log ( ''tijX;^^ ) - ^ t + oig% (15) 

Next we consider the pure scalar part of the action. Introducing the notation 

^Px = Ho{x) + iTrTTr{x), (16) 

it reads: 

Sh — E {^0(2;)^ + TTr{x)Trr{x) + X{Ho{x)^ + TTr{x)Trr{x) — i f 

X 

3 

-2ks ^ {Ho{x + fia^)Ha{x) - 7rr{x + //a^)7r,.(a::)) 

{x + lai)T:r{x))] . (17) 
Assuming that the i?o field has a non-zero vacuum expectation value we write: 

Ho{x)^H{x)+v. (18) 

Moreover we introduce the notations 
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and the continuum fields 



5o = 6A^^, A, = ^ (19) 



/fe(x)=f^)'F(a;), <(a;)=(^)\'-(a.), = f ^ ) v. (20) 
\asatj \asatj \aaat,' 

Using these we find the scalar part of the action using continuum variables to be 
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Sh = alat ^ \ \ Y,{ViH,{x)ViH,{x) + Vi<(a;)Vi<(a;)) 



where 



and 



2 , 

+ ^{ViHc{x)V4Hc{x) + V4<(X)V4<(X)) 

+ml{mx) + v^f + <(x)<(a.)) + ^m{x) + v,f + <(a;)<(a;))2} , (21) 
alml='-^-Q-2^l (22) 



V,f{x) = + (23) 



is the lattice derivative. 

Putting Vc = Q above corresponds to the symmetric phase, in this case mg > 0. Determining Vc from the non-trivial 
minimum of the scalar potential one gets 



2 



6m: 



2 



^2 _ 



^ for ml < 0. (24) 



Introducing mj^ q = — 2mQ we obtain finally: 



Sh = alat ^ ^ i ^(V,i/c(xOV,i/c(x) + Vi<(x)Vi<(x)) 

xeA li=l 



7^, 



+^{ViHc{x)ViHo{x) + V4<(a;)V4<(a;)) 
+™2^,oif,(x)2 + g(ife(x)2 + <(x)<(x))2 + ^JJ,,o(a;)(i?,%(x) + <(x)<(x))} . (25) 
Now we consider the gauge-scalar interaction: 

xeA l i=l J 

Introducing continuum variables we obtain 

Si = alat V I V(a°(x) - 1) {vl + 2^;eifc(a;) + H,{x)H,{x) + <(x)<(x)) 

xeA U=i ^ 

-^ViHcix) - — {Hc{x)ViHc{x) +Trlix)ViTr^^{x)) 

+7f K(a;) - 1) {-\ {vl + 2vcH,{x) + Hc{x)H,{x) + <(xX(x)) - ^ViHdx)- 
i; \ ttj at 

— (ire(x)V4i?c(^) +<(a;)V4<(x)) 
do 



+ -J2<{x) (e„t<(a;)Vi7r*(a;) + 7rl{x)ViH,{x) - {H,{x) + Vc)Vi<(a;)) 
+^jlal{x) (e,,t<(a;)V47r*(x) + <(a;)V4ffc(a;) - (i?c(a;) + ^;c)V4<(x)) 1 . (27) 



In perturbation theory the gauge has to be fixed. We use as the gauge fixing function 

A al{x) -a[(a;-m,) -f^{al{x) - al{x - 4at)) aVcg^TT^c 

1=1 ' * 

which is a lattice version of the well known continuum gauge fixing function. In eq. ( p8|) a is the gauge parameter. 
This choice ensures that the mixed second order term in A^^{x) and tt'^{x) will drop out from the sum of the gauge-scalar 
interaction and the gauge fixing parts of the action. We obtain 

= E fr{^)fr{^)^ (29) 



^ r 1 

Sfp = a^at X! X! | ~2 - Cr{x + iai))al{x)(cr{x) + Cr{x + iai)) 

xeh 1=1 ^ 

+erst(cr(a;) + Cr(a; + iai))al{x){ct{x) - Ct{x + iaj))| 



„2 



+ 7T% {(cr(a;) - Crix + 4:at))al{x){crix) + Crix + lat)) 
+£rst{crix) + Cr{x + 4:at))al{x){ct{x) - Ct{x + 4a4))} 
-a^-^ {cr(a;)cr(a;)(-ffc(a;) + Uc) + erstCr{x)cs{x)'Kl{x)] > . (30) 



The final form of the continuum notation action reads: 

Scant — Spf^ + S'pi'''^ + Sm + Sh + Si + Sgf + SpP, (31) 

where the individual terms are given in eqs. (|lj, ^ |l5|, |2^, 29, ^0|). The vertices of perturbation theory may be 
easily obtained from eq. (^). As usual in lattice perturbation theory we have new vertices proportional to for 
n > 2. 

The formulae to compute the Fourier transforms are as follows. For the gauge field: 

«fe,/x = alat exp ( -i{k, x) - ^fc^ ) a^(x), (32) 



where {k,x) = 27r(^ + ■ •• + ^), i? = 0, 1, 2, 3 and 



moreover the integers i^^, (Z^) take values from 0, 1, . . . , — 1. 
The inverse relation is: 

= T-TTl E '^^P ( *(^' ^) + ^^A. ) ■ (34) 

L1L2L3L4 aia2a3a4 V 2 ' 

For a lattice infinite in all directions we have 



1 TT- /■ ' 

(a;) = n y / ^^Pp exp (ip • a; + ip^a^) af^^, (35) 



and p- X ^ Y.t,=iPvXy. 

For scalar fields we have similar formulae, however, the second terms are missing in the exponents of eqs. (32 3^j35| ). 

Our aim is to perform perturbative calculations. The first step is to write down the propagators of the fields from 
parts quadratic in the fields of the action. We want to determine the tree-level propagators, which are zeroth order 
in g and A. Since 7/3 and 7^ do depend on the couplings, we use ^ in the propagators as their tree-level values. The 
remaining correction terms from 7^ and 7^ are quadratic in the fields and give two particle vertices similarly to the 
measure term in the action in the isotropic case. These will be absorbed by the kinetic parts of the propagators (see 
later in eqs. (||)-(|7|)). 

The inverse tree-level propagators in momentum space have the following forms. 
For the Higgs boson 
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(36) 



for the Goldstone bosons 



(37) 



for the gauge boson 



i=l 



P^Pu- 



1 + a 



(38) 



for the ghost 



Afp,o(p) ^ ^^Pi ~ am^Q, 



(39) 



where 



2 . asPi ^ 2 . atPi 
= — sm— -, p4 = — sm— -. 



(40) 



The masses have the fohowing expressions in terms of other parameters: 



1 - 2A 



-C - 6 - 2^^ 



, '"wo 



2AC/3 4 ' 



(41) 



B. Critical hopping parameter and anisotropy parameters 



The main goal of the paper is to perform a 0{g^, A) analysis of the theory defined by eq. (|l]). This means first the 
determination of the mass-counterterms. One wants to tune the bare parameters in a way to ensure that the one-loop 
renormalized masses are finite in the continuum limit (however, their values in lattice units do vanish, agniren — 
for — > 0, a* — > 0, ^ — fixed). At the same time the vacuum expectation value of the scalar field will be also 
zero in lattice units (a^Uc = 0), i.e. we are at the phase transition point between the spontaneously broken Higgs 
phase and the SU(2) symmetric phase. The condition is fulfilled by an appropriate choice of the hopping parameter 
(critical hopping parameter). The ratios of the couplings (7^ and 7^) are still free parameters and can be fixed by 
two additional conditions. We demand rotational (Lorenz) invariance for the scalar and gauge boson propagators on 
the one-loop level. This ensures that the propagators with one-loop corrections have the same form in the z- and t- 
directions. Clearly, arbitrary couplings for different directions in eq. (|l|) would not lead to such rotationally invariant 
two-point functions. 









FIG. 1. Higgs boson self energy graphs. Solid lines stand for Higgs, dashed lines for Goldstone, curly lines for vector bosons 
and wavy lines denote ghosts. 
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The most straightforward method to determine the transition point is the use of the effective potential. The condition 
(PVeffi^ = 0)/d^^ — gives a simple, gauge invariant expression for the value of the critical hopping parameter, 
which is exact in the continuum limit. The relevant formulae are given in Landau gauge in [p3| . Here we present the 
general i?^ effective potential: 



K//(<I>) = ^$2 + A,<1>^ 



24A, 



■ log(A: 



16Ac — ag'^ 
+61og(P + 52$V4) 



5 P^ag^ 
2 I A:2 + 52$2/4 



(42) 



where 



(27r) 



(43) 



Alternatively, one may calculate the one-loop corrections to the masses and require that the renormalized masses be 
zero in lattice units in the limit of zero lattice spacing and fixed ^ = ag/at, as explained above. 

First we consider the corrections arising from the two-point interaction vertices. In addition to these there are the 
one-loop corrections, which we evaluate later on. Including the two-point interaction vertex corrections the momentum 
squared sums in eqs. (36), (Is^), (^) modify to 



i=l 



Pt ' Z^Pz + -72P4. 



(44) 



for the Higgs, Goldstone and ghost propagators. The gauge boson inverse propagator becomes more complicated: 



~^t,,{pr' 



9 9 



PzPj 1 + ^ 



(45) 



21 



9sgt 7k 

ag^ ^2 



(46) 



9993 24 

"c5t 7k , ~2 -2 9t 7k 



(47) 



Let us now consider the self-energy corrections to the Higgs mass. The relevant diagrams are shown in figure 1. 
Evaluating all graphs we obtain at zero Higgs four-momentum, independent of the gauge choosen (i.e. independent 
on a): 



a2(mg)2 ^ alml -(2go + ^g' ] M^,0), 



(48) 



where we used the notation: 



J„{^,mas 



4-2n 4 rTT/ap 



n 



dkn 



(27r)4 (m2 + fc2)n 



(49) 



Inserting 



(^Irnjj = -2 



1 - 2A 



- 6 - 2^ 



(50) 



for the one-loop corrected bare mass and using the notation A = K'^go/{6^) = Ak^Xc/S, together with a2(m|^)2 = 0, we 
get solving perturbatively for k: 



1 



2(3 + ^2) (3 + ^2)2 



6CJi(e,0) 



9gJi(e,o) 2 

16(3 + ^2)2^ ■ 



(51) 



This result coincides with the d?Veff{^ ~ 0)/d^^ = condition of eq. (12) of For the readers' convenience we 
plot Ji(^, 0) of eq. ( |49| ) in figure 1 as a function of 1/^. For the special case of symmetric lattice spacings, ^ = 1, our 
quantum corrections to the critical hopping parameter reproduce the known result of the isotropic SU(2)-Higgs model 
( ||23,|2^). An 8-term Chebishev polynomial approximation with 6 • 10~^ accuracy to the function reads: 



Ji(^, 0) 0.2276734 - 0.000175561/^ - 0.1452559/^^ 
-0.03593908/C^ + 0.3487585/^"* - 0.4128226/^^ 
+0.2187872/^6 - 0.04609285/^^. 



(52) 



It is instructive to check that the same result is obtained starting from the symmetric phase perturbation theory, 
where some graphs are absent and one is lead to 



= — a^TOn 



9a 



1 Ji(e,0). 



(53) 



0.22 - 



0.2 



0.18 



0.16 




0.4 0.6 

FIG. 2. The lattice integral Ji{£,, 0) (see text) on asymmetric lattices as a function of 1/^. 



Let us now consider the self-energy corrections to the gauge boson mass. The relevant diagrams are shown in figure 
3. The inverse propagator (eqs. (^5|)-(|4^)) at zero momentum has a specific structure, namely 

A?^..,(0)-'-rt,m2,4' 
9 

A^.44(0)-^ = <5«^m2,^^. (54) 

One therefore has to determine both the diagonal space-space and the time-time components in order to check 
consistency. Since the bare mass squared turns out to be 0{g^), we may safely put gl/g^ = g\lg^ = 7k/^^ = 1 ™ 
(^). Finally we obtain, after imposing zero rcnormalized lattice unit mass squared: 

2 2 2 ( , ^ ™w 
c " ' 

Inserting 



a>^=.9'( 2 + 2;^Mi(^-0)- (55) 



2 2 ira\ 

m^Y^mfj- — , (56) 



450' 

we get back eq. (^ij) consistently. Again we have checked that (^^ holds in all gauges 
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Next we discuss the anisotropy parameters 7/3 and 7^. Following Karsch and Stamatescu we determine them 
from the requirement of rotational invariance in the continuum limit as, at at fixed ^ — as/at- In particular we 
consider the physical particle propagators, which receive quantum corrections 

AnAPr' - A«(p)-i + S^ab), Aw,i,M^b)"' = ^w.^^iPr' + ^w,i,,M, (57) 
where Ah{p)~^ and Aw,fii/ are the tree-level propagators corrected with the two-point vertices. Ah{p)^^ is given by 

A^(p)-i=m^,o + ^p? + JpI (58) 

while Aw,p.u is given by (||)"@. 




The corrections to the anisotropics in the kinetic parts of eqs. (^) should be cancelled by the kinetic parts of the 
self-energies. For the Higgs boson this can be achieved by requiring 



1 



2 



dpi 



p=0 



e 



dpi 



(59) 



where z = 1, 2, 3. The graphs contributing are the momentum dependent ones of figure 1. 
For the gauge boson there are several possibilities. As a simple one we choose 



1 d^'Ew.iAp) 
^2 dp} 



p=0 



l^^^w^M{p) 



dpi 



p=0 



(60) 



where i j — 1,2, 3. This is easily calculated, since only the Sij term of A^^j(p)^^ contributes on the left hand side. 
Not all the self energy graphs contribute, but only those graphs of figure 3, which depend on the momentum. 
Our results for infinite lattices are 



(61) 



/■"" /■"" /■"" ,4 rr/ 1 \ 3 + cos((7iai) cos((73a3) - 003(9404) 

'"^^^ ^ J-^ '-^ '-^ 'mnv'e) — i i 

- ^2 ^ ^2 (^^ ^ ^ j sin^((7iai) + ^ cos^(^yi) (003(9303) - 003(9404)) 
y\ / sin^(g3a3) - sin^(g4a4) ^ 003(9303) - 003(9404) 

M „2 o 



8sin^(9iai) 



2co32(aiiii 



sin^(9iai) /sin^(93a3) - sin^(94a4) ^ 

+ I ^fEl 003(93.3) + 003(9404)^ 



(62) 
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4 7_^y_^y_^y_^^^E'?, 



16 



(^^sin^(g4a4) - sm^{qiai)) I , 



(63) 



where the sums are over fi = 1, ... ,4. The above expressions are easily seen to be finite and independent of Os and 
at- We have also checked that they are gauge independent. The dependence on ^ is plotted in figure 4. 
A 6-terni Chebishev polynomial approximation with 2 • 10~^ accuracy to the functions reads: 



c/3(0 = -0.1687249 + 0.124013/C + 0.08608489/C^ 
-0.04715295/^^ - 0.0002526438/^* + 0.006038775/^^ 



(64) 



c^{^) = -0.05691582 - 0.0001275536/^+0.07582766/^^ 
-0.003112956/e^ - 0.0265274/^'' + 0.01085953/^^ 



We also have to equate Avy,i,i3 and Avk,i,i4 

1 , 9. 



dpidpz 



p=0 



p=0 



(65) 



(66) 



This is a non-trivial constraint, which our previous expressions do satisfy. 



-0.05 



-0.1 



-0.15 




0.4 0.6 
FIG. 4. C/3(5) and Ck(5) as functions of 1/^. 



There are several important features of the anysotropy parameter result, which should be mentioned. 
a. Masses in the propagators: A consistent perturbative procedure on the lattice determines the bare parameters, 
for which the renormalized masses vanish, cf. eq. (|5l|). With these bare couplings other quantities, e.g. asymmetry 



parameters, are determined. However, using the one-loop renormalized masses (agm^ 



agrn 



w 



0) in the propagators 



instead of the bare ones leads to changes in the results, which are higher order in g^ and A. Therefore, all our results 
are given by the integrals with renormalized masses. 

h. g^ and A corrections: In figure 4 we have given only C/3(^) and Ck(^). As shown by eq. ( ^l|) the functions 6/3(0 
and 6k(0 vanish, thus there are no corrections of 0{X) to the anisotropy parameters. It is easy to understand this 
result qualitatively, since only graphs with two or more scalar self-interaction vertices have non-trivial dependence on 
the external momentum. This feature is connected with the well-known fact that the theory does not have any 
wave function correction in first order in the scalar self-coupling. It is worth mentioning that there is only one type 
of two-loop graph (the setting-sun) which should be combined with the one-loop graphs, in order to obtain the whole 
O(A^) correction. 

c. Pure gauge theory: A number of graphs of figure 3 (namely those containing only vector boson and ghost lines) 
are identical to those of the pure gauge theory. Evaluating the momentum dependent ones from these diagrams, one 
reproduces the result of ref. ||2l[ (the function C/3(0 of the present paper corresponds to Ct{C) — Cct(C) of ref. |^). The 
most important contribution comes from the self-energy graph with gauge boson four- coupling. Inclusion of the scalar 
particles gives only small changes. The relative difference between the cp{S^) functions for the pure SU(2) theory and 
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for the SU(2)-Higgs model is typically a few %. 

d. Quantum corrections to the hopping parameter: The contributions to the hopping parameter come from the 
momentum dependent graphs of figure 1. This correction has the same sign and order of magnitude than that of 
the gauge anisotropy parameter; however it is somewhat smaller. It is possible to combine the anisotropics c'^^iO = 
CfsiO ~ Ck(0- Foi' this choice in the gauge sector and with 7^ = ^ the rotational invariance can be restored on the 
one-loop level, choosing the appropriate value for the lattice spacing asymmetry as /at- Thus, the masses in both 
directions will be the same. However, the obtained lattice spacing asymmetry will then slightly differ from the original 
e One gets as/ at = C(l - .9'ck(C)/2) + 0{g\ A^). 

e. For later use we specify: cp{4:) = -0.13308, c«(4) = -0.052353, thus c^(4) = -0.080727, 7^(4) = 3.9193, 
as/at = 4.05235. 

/. Asymmetry parameters away from the critical line: Following the procedure outlined above one may determine the 
asymmetry parameters away from the critical line. In this case tree-level masses are nonvanishing and are in fact 0(1). 
Therefore one has to keep them in the propagator denominators. Thus the final results become more complicated. We 
do not reproduce the formulae here, only note that numerically the results are very close to the previous case. Thus 
the asymmetries determined near the critical line are universally applicable. 

g. Finite lattice results: The above formulae are valid for infinite lattice sizes, however, replacing the lattice integrals 
with the appropriate lattice sums, one gets results valid for finite lattices. 



C. Perturbative study of the continuum limit of the finite temperature theory and optimal choice of the 

parameter ^ 

The approach to the continuum limit of the finite temperature theory may be studied in the approximation of one- 
loop perturbation theory. The relevant physical quantities we study are the ratio of the critical temperature (Tc) and 
the Higgs mass and the ratio of the Higgs and vector boson masses. To calculate them in perturbation theory we first 
determine the bare Higgs mass parameter using the analogue of eq. (|4^) for a lattice with finite extension (Lt) in the 
^-direction, i.e. at finite temperature T = l/{Ltat), by imposing the condition a1{m^)'^ = 0. This choice corresponds 
to the lowest point of the metastability region with = l/{Ltat), i.e. when the derivative of the effective potential 
at zero field first becomes negative. Using the same bare coupling parameters in the action we next determine the 
physical Higgs and vector boson masses on a T=0 lattice (i.e. using a lattice with equal (infinite) physical dimensions 
in space and time directions). 

More precisely, the bare quantity a^m^ is determined from Eq. ( ^8|) with a^(m||)^ — 0, replacing however Ji(^, 0) 
with JT{Lt,£,,0), where 

JT{Lt,£,,mas) = —^\{ dkp (67) 



and in the denominator is given by 



ki — — sm . (68) 



at Lt 

(It is straightforward to write down the finite lattice version of eq. (p7|), too.) The T = renormalizcd Higgs mass 



(a^m^l) is then determined from the unmodified eq. (48) using the already known value of the bare parameter a^m^ 
and the infinite volume T = integral Ji(^,0). Using Tc — \/{atLt) = ^/{osLt) we finally obtain the simple formula 
for a given Lt: 



Tc/m^ = 4- : -. (69) 

Lt[{2go + 9gy2){MLt,tO)-Ji{^,0)P 



In the same approximation {m^/m^r)^ equals to the tree level value 4:go/{3g'^) 
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FIG. 5. Simulation time necessary to reach 0.1% precision determination ofTc/m^ (normalized to the ^ = 1 point) versus 

The result eq. (|6^) refers to infinitely large lattices (i.e. infinite in both the space-like and time-like directions for 
the T = case and infinite in only the space-like direction for the T ^ case.) The continuum limit (realized as 
Lt oo) is well defined. In lattice simulations, however, we always have finite lattices. We have to choose minimal 
lattice volumes large enough to ensure a reasonable precision. This choice of course does depend on i^, therefore we may 
also look for the optimal choice of ^ ensuring a reasonable precision (say 0.1%) of the physical mass determinations 
using the smallest possible lattices or shortest simulation times. This problem may be studied in lattice perturbation 
theory. 

To obtain the optimal choice of ^ we first determine the Lt ^ oo (i.e. the continuum) limit value of Tc/m^ as a 
function of ^ using eq. (|9|). We obtain that - as expected - the limit of Tc/m^ does not depend on ^ within errors. 

Next we take into account that in practice we simulate on lattices with finite extensions. In order to fit in the 
relevant modes we have to deal with a given physical volume: 

W = Ltat{L,Catf = ^{L./Lt ■ i/Tf. (70) 
Thus the number of the lattice points (which determines the memory required) is expressed as 

LtLl = VphysT^Lt/e- (71) 

To get a correct estimate of the simulation time we have to take into account the autocorrelation times as well. Since 
these are proportional to the squares of the correlation lengths for a local updating algorithm (see [^), i.e. to {Lt/£,)^, 
the time necessary for simulation on a given physical volume and temperature will be proportional to 

Vp^,y,T^L^Je■ (72) 

Next we choose a lattice extension in temporal direction Lt so that by eq. ( |69| ) we obtain an approximation of the 
previously determined continuum limit Tc/m^ value to a given (say 0.1%) precision. Tc/m^ as determined from 
eq. ( |69| ) as a function of Lt approaches the limiting value from below for large Lt for all ^ values. However, for 
^ > 2 it decreases for increasing, small Lt values. Thus specific small Lt values may better approximate the limiting 
value of Tc/m^ than larger intermediate values. It is clear that this is an accidental agreement only, therefore in our 
considerations we have determined the smallest Lt value giving Tc/m^ with the required precision, which does not 
deteriorate for larger Lt. 

More precisely we compare the true continuum limit of Tc/m^ with an approximate value obtained from an ex- 
trapolation to Lt = oo of the Tc/m^ values determined from four subsequent Lt values. We choose Lt to be the 
minimal Lt, which (together with the 3 larger Lt values) already gives the required precision. Having determined Lt 
we calculate the corresponding simulation time for finite lattice size using eq. (|72|). Figure 5 shows the simulation 
time normalized to the ^ = 1 value as a function of 1/^ for 0.1% precision in Tc/mM. The normalized simulation time 
as a function of ^ has a broad minimum near ^ = 2. The number of lattice points (|7l] ) (normalized to the ^ = 1 value) 
is quite a similar function of ^ with a broad minimum near ^ = 2. In our numerical simulations we have chooscn 
^ ~ 4, which is a good choice both from the point of view of simulation time and filing in the relevant modes into a 
practically accessible lattice. 
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III. NON-PERTURBATIVE ANALYSIS OF THE ANISOTROPIES 



This section of the paper deals with our non-perturbative determination of the anisotropy parameters by means of 
numerical simulations. Besides a mere confirmation of the one- loop calculations in the previous part, it could give 
estimates of possible corrections, which go beyond perturbation theory. This is an important step towards future 
studies of the finite temperature electroweak phase transition in the framework of the four-dimensional SU(2)-Higgs 
model on anisotropic lattices. Namely, if the deviation from the perturbative results turns out to be so small that 
its influence on expectation values in a numerical simulation is negligible within their typical statistical errors, the 
one-loop perturbative anisotropics 7/3, 7^ and ^ can be used without any further (non-perturbative) fine-tuning. At 
first sight this may not seem very surprising, because the zero-temperature theory is weakly coupled {g^ ~ 0.5). But 
owing to the fact that the corrections in the parameter A — entering only at two-loop level — whose size essentially 
determines the value of the Higgs boson mass, are not exactly known, such an investigation is necessary, particularly 
in view of Higgs masses around 80 GeV or larger, which is the physically allowed region determined by the LEP 
experiments. 

As already discussed above, the tree-level values of the anisotropics receive quantum corrections, which in general 
have to be determined non-perturbatively. A physically motivated idea for their estimation is to impose the restoration 
of the space-time interchange symmetry as a remnant of Lorentz invariance after discretization of the continuum theory. 
In practice this is to be realized by the requirement that Higgs and gauge boson correlation lengths in physical units 
should be equal in space- and time-like directions. Furthermore, we include into the analysis the length scale of the 
static potential derived from space-time and space-space Wilson loops. 

The following subsections describe our numerical studies in more detail. After some brief remarks on the simulation 
techniques and parameters used, we present the results on the physical observables under consideration and propose, 
how they can serve to extract the coupling and lattice spacing anisotropics non-perturbatively. Finally, the values 
obtained in this way are confronted with perturbation theory. 



A. Monte Carlo simulation and its parameters 



In our Monte Carlo simulations we apply an optimized combination of heatbath and overrelaxation algorithms, 
which has been extensively discussed for the isotropic model in refs. p^ , p^ , p8[ | , and their implication carries over 
straightforwardly to an anisotropic lattice. The action (|l|) is easily arranged to S\U,ip\ = X^ajGA'^^" ^'^'^ lattice 
action per point 

5. = 6/3Pp,, + i?, + AQ, - 8Ki^,, (73) 
consists of the length variables of the Higgs field 

i?, EE iTr (^+^,) = Q, = {pl-l)\ (74) 

of the weighted sum of the plaquette contributions Ux-.iiv = Ux^f^Ux+fi.ijU^^- iJ^x,v ly™S the space-space and the 
space-time planes 



Pp,s,.^l (l-iTrf/,;^,), Pp,t,.^^ (l-iTrC/.;;..), (76) 

and of the weighted sum of the space- and time-like components of the (^-link operator L^p-x^ = ^Tr (^(p^^fj^l/x.^ii^x^' 

1/3 \ 

Lm^X = T ( Lip_s,X + ^K,Lip^t,X 1 (77) 

4 V7k J 



For ^ = 7/3 = 7k = 1 this action simplifies to its well known form on isotropic lattices. Eqs. (|7^ ) - (^8|) already cover 
most of the observables, whose expectation values are calculated by numerical simulations. 

The updating scheme per sweep, a sequence of one Ux,^- and one (^2;~heatbath step, succeeded by one Ux,^- and 
three <pa;-overrelaxation steps, has been taken over from refs. 0,p8l. There it was observed that the inclusion of the 
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overrelaxation algorithms reduced the autocorrelation times substantially, in particular for the operators and 
Lip^ whose expectation values show the largest autocorrelations. 

As pointed out in the introduction, the anisotropic version of the SU(2)-Higgs model is believed to provide quantita- 
tive insights into the electroweak phase transition at large Higgs boson masses of tuh > 80 GeV, at which the typical 
excitations with small masses (i.e. large correlation lengths) would demand very large isotropic lattices exceeding any 
presently accessible computer resources. In principle a rough resolution in the spatial directions by moderate lattices 
combined with accordingly large lattice anisotropics ^ could handle this situation. However, for T > a small temporal 
extension Lt sets the (very large) temperature scale through T — 1/atLt, and hence it is more sensible to ensure a 
large enough lattice cutoff by employing ^ Lt, thus in our numerical work we take 

C-4, (79) 

which is also strongly motivated by the result of subsection III.C. Since this makes the correlation lengths in time 
direction smaller than in space directions, it seems to be reasonable to fulfill Lt ~ S,Lz in order to restore the symmetry 
of the physical extensions and to enable a precise mass determination. We consider two lattices of sizes 8^ x 12 x 48 
and 8^ X 16 X 64, where the spatial correlation lengths correspond to few lattice units and the finite-volume effects are 
expected to be small. 

The T ~ simulations are generically intended to fix the physical parameters, i.e. renormalized couplings and 
masses. Consequently, the lattice parameters in this study are chosen to reach the interesting region of mn — 80 GeV 
or a Higgs to gauge boson mass ratio of 

Rhw = — ^ 1 (80) 

with the experimental input mw — 80 GeV setting the overall physical scale. This is (at least approximately) achieved 
by the values f3 = 8.0 and A = 0.000178. The scalar hopping parameter, which has to comply with the condition 
that the T > system is at a phase transition point for a certain temporal lattice extension, is calculated from the 
discretized version of eq. (|l|)[|. Referring to if = 4 this amomits to k = 0.10662. The non-perturbative corrections 
usually tend to decrease the tree-level mass ratio 



Rhw,o = — ^ = \ — 1^ • (81) 
mw,o V K 

Our strategy for the determination of the coupling anisotropics is as follows. In the numerical simulation we have to 
find those couplings of eq. (|l|), for which the space-time symmetry is restored. Therefore, we fix one of the coupling 
anisotropics to its tree-level value, ignoring its quantum corrections, and tune the other one to produce identical ratios 
of (decay) masses in space- and time-like directions for a set of two or more (particle) channels. The mass ratios 
determine the actual lattice anisotropy, which will then slightly differ from the original ^ of (|79|). In this spirit we 
choose three pairs of coupling anisotropics, denoted as 'tree', 'low' and 'perturbative'. 



t : 


7k 


= 4.0, 


7/3 


= 4.0 


1 : 


7k 


= 4.0, 


7/3 


= 3.8 


P : 


7k 


= 4.0, 


7/3 


= 3.919, 



(82) 

and calculate the corresponding lattice spacing anisotropics from different physical quantities as described compre- 
hensively in the subsequent subsections. Assuming that they depend linearly on 7^ in this small interval, we can 

interpolate to a matching point (j^^^K S,^"^^"^) , at which all ^-values coincide within errors. These estimates are quoted 
as our non-perturbative results. 

All numerical simulations have been done independently on the APE-Quadrics computers at DESY-IfH in Zeuthen, 
Germany, and — to a smaller extent — on the CRAY Y-MP8 and T90 of HLRZ in Jiilich, Germany, which offer 64-bit 
floating point precision. In contrast to some quantities, e.g. the critical hopping parameter in T > simulations, the 
32-bit arithmetics of the APE-Quadrics is sufficient for the calculation of all T = quantities, especially for particle 
masses and the static potential. 

B. Correlation functions and masses 

We now turn to the determination of the Higgs and gauge boson masses. As in refs. |l|,|l|, they were obtained 
from suitable correlation functions of gauge invariant, local operators integrated over time (space) slices. Those are 
Rx and L^-^fi for the Higgs mass, and the composite link fields 



^ The knowledge of the more accurate, non-perturbative value of the critical hopping parameter, which has to be determined 
numerically, is not relevant here. 
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x:rk 



= -Tr (^Tra^^j^Ux.kOix^ , Tr- Pauli matrices , r,k = 1,2,3 (83) 



for the gauge (M^-boson) mass. 

The connected correlation functions Tq of these operators have been measured in the time-like and in one space-like 
direction. For the Higgs mass mn the functions To{t) and ro(z) were calculated from t~ and z-slice averages of Rx 
and the weighted (p-link L^^^^ of eq. (^7|). Since these functions can not be regarded as uncorrelated, we have averaged 
them — after an appropriate normalization of the correlations at distance zero — before performing the mass fits. 
The same prescription holds for the gauge boson mass mw, but with two major differences: firstly, the t- and z-slice 
correlation functions of Wx;rk have been measured separately for all combinations of r and k, and secondly, in place 
of fc = 3 in (|3|) actually we have to take fc = 4 for the correlations in z-direction (i.e. all directions in Wx-rk are 
orthogonal to the direction of propagation). Again the individual correlation functions are averaged to one function 
per direction as in the Higgs channel. 

As lowest energies the particle masses are extracted from one-exponential least squares fits to shapes of the form 

ro(^) = ^[e-"^ + e-'"(^-^)] +C, e^O,l,...,^, Le{Lt,L,^ (84) 

with m G {atmH,t, o-tfnw,t} or m S {asmn^s, cis'mw,s}, respectively. The constant terms in the vector channel are 
highly suppressed so that a two-parameter fit is mostly sufficient. Each full data sample has been divided into 
subsamples, and the statistical errors on the masses originate from jackknife analyses. All simulation parameters and 
lattice sizes are collected in table |. 

Our fitting procedure consists of correlated fits, sometimes with eigenvalue smoothing, and simple uncorrelated fits. 
For the former we use the Michael-McKerrel method j3^, whose features and application in the SU(2)-Higgs model 
have been sketched in ref. [Q. Its main purpose is to select the most reasonable fit interval in data sets, which are 
strongly correlated in the fitted direction. Uncorrelated fits, which ignore these correlations, are often plagued with 
very small values of P^r degree of freedom (dof) for nearly all fit intervals in question, whereas in correlated fits 
the emergence of x^/dof ~ 1 for some fit intervals represents a safe criterion to select reasonable fit intervals. This 
also works well for data sets of lower statistics, if the smallest eigenvalues of the correlation matrix are smeared via 
replacing them by their average. All resulting mass estimates in lattice units are shown in tables ^ and III. We chose 



the largest fit interval with a reasonable x^/dof from the correlated fit and the results of the uncorrelated fit along 
this interval as the final fit parameters. Both fits were always consistent within errors, and other fit intervals with 
comparable or even lower x^/dof did not cause any significant changes. 



index 


lattice 


7-^ 


7/3 


correlation functions 


Wilson loops 


sweeps 


subsamp. 


subsamp. 


indep. sweeps 


tl 
11 
pl 


8^ 
8^ 
8^ 


X 12 X 48 
X 12 X 48 
X 12 X 48 


4.0 
4.0 
4.0 


4.0 
3.8 
3.919 


100000 
100000 
576000 


50 
50 
192 


50 
50 


100 
100 


t2 


8^ 


X 16 X 64 


4.0 


4.0 


192000 


64 


64 


150 


12 


8^ 


X 16 X 64 


4.0 


3.8 


192000 


64 


64 


150 


p2 


8^ 


X 16 X 64 


4.0 


3.919 


704000 


256 


128 


150 



TABLE I. Summary of the numerical simulations for the mass and static potential computations. The other parameters are 
f3 = 8.0, A = 0.000178, and k = 0.10662. 
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As emphasized above, the space-time symmetry restoration, which implicitly establishes becomes apparent in 

equal physical correlation lengths as = at(, of the theory. Thus we introduce anisotropy parameters in the Higgs and 
vector channels by calculating the ratios 

_ asmH,s ^ _ asmw.s ,oc^ 

4h = , i.w = [ob) 

atmH,t atfnw.t 

within the jackknife samples of the space- and time-like masses. These are displayed again in tables and Due 
to the compatibility of the results from the two lattices one concludes that the finite-size effects are quite small. 



quantity 




tl 


11 


pl 


atmH,t 
asmH,s 


4 - 18 : 0.1408(22) 
1 - 6 : 0.5635(31) 


4 - 22 : 0.1370(27) 
1 - 6 : 0.5611(62) 


4 - 24 : 0.1387(15) 
1 - 6 : 0.5603(30) 






4.002(67) 


4.097(86) 


4.041(55) 


a,stnw,s 


8 - 24 : 0.1523(13) 
1 - 6 : 0.6225(29) 


8 - 22 : 0.1538(13) 
2 - 6 : 0.6066(40) 


8 - 24 : 0.1554(25) 
2 - 6 : 0.6307(22) 






4.091(30) 


3.945(32) 


4.059(44) 


RHW,t 

Rhw,s 




0.925(15) 
0.905(6) 


0.891(20) 
0.925(12) 


0.892(18) 
0.888(7) 



TABLE II. Fit intervals, Higgs and gauge boson masses in time- and space-like directions, and the resulting lattice spacing 
anisotropies for the smaller lattice. 



quantity 




t2 


12 


p2 


atmH,t 
asmH,s 


4 - 32 : 0.1408(22) 
1 - 8 : 0.5590(42) 


4 - 32 : 0.1370(27) 
1 - 7 : 0.5586(40) 


4- 32 : 0.1378(11) 
1 - 8 : 0.5550(40) 






3.969(73) 


4.078(80) 


4.027(36) 


a,tm.w,t 


8 - 32 : 0.1499(31) 
1 - 8 : 0.6318(40) 


8 - 30 : 0.1599(42) 
3 - 8 : 0.607(11) 


6 - 32 : 0.1525(15) 
2 - 8 ; 0.6133(27) 






4.23(10) 


3.80(13) 


4.021(48) 


RHW,t 

Rhw,s 




0.940(24) 
0.885(8) 


0.857(26) 
0.921(20) 


0.904(11) 
0.905(5) 



TABLE III. The same quantities as in table O for the larger lattice. 
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C. Wilson loops and static potentials 



Another approach to the ^-determination is based on the static potential, which has the physical interpretation as 
the energy of an external pair of static charges brought into the system. To this end we have measured rectangular 
on-axis Wilson loops Wij{Ri, Rj) of extensions 1 < -Ri < Li/2 and 1 < Rj < Lj/2, lying in space-time and space- 
space planes. The gauge configuration was transformed to temporal gauge for space-time and to A^^^x) = gauge 
for space-space Wilson loops, and every loop with two sides in t- or z-direction, respectively, was included in the 
statistics. 

As a generalization of the isotropic lattice case we distinguish between static potentials 

lim In W,,{R.„R,) (86) 

in space-like {ij — st, ss) and time-like [ij — ts) directions, according to the Rj oo extrapolation in the second 
argument of Wij, which is supposed to be done first. The shape of the potential, which is governed by a massive W- 
boson exchange is known to be Yukawa-like, and calculating along the lines of refs. [ p9|j30| lowest order (tree-level) 
lattice perturbation theory yields 

with lattice momenta fc„ — 2a~^ sin (a„/c„/2), n = 1, . . . , 4. In the continuum limit this expression reflects the usual 
screening behaviour, i.e. modulo a constant, 

independent of i and j. After substituting p„ — Unkn with p„ = 27r^„/L„ and /„ = 0, 1, . . . , L„ — 1 on a finite lattice, 
one obtains from eq. ( ^7|) 



35^ ^ 

where mv,ij = aimw,o and 



IDTT 



T t o ^ - ^ cos {R,pi) 

I^J{mv.^J,R^) = —— }^ -^—^ -^^ . ^ ^ r, (90) 

where k and I are different from each other and from i and j. 

Since = gji + 0((?^), the simulation results for Vij are fitted with the ansatz 

a^V,J{R,) = e-™^-^' + a, + D,jG,,{mv,^,Ri), (91) 

Ri 

where is a term correcting for finite-lattice (size and spacing) artefacts, and Aij, my ij, Cij, Dij are the parameters 
to be fitted. Gy reads: 

Gy (my.„- ,R,)^^e - J,^- {mv,^J ,Rr). (92) 

Hi 

By definition the " global" renormalized coupling is obtained by identifying the coefficient of the contribution relevant 
at short distances, 

= (93) 

Note that rnv,ij/ai and also g'^ as determined from Wilson loops with different indices have to be independent of the 
indices for properly choosen coupling anisotropics. 

In a first step of the analysis we performed multi-exponential fits Wij{Ri, Rj) — X]^=o ^ order to get 
the potential for fixed Ri as the ground state energy Vq from the large Rj asymptotics of the Wilson loops in (p6|). 
Starting at distances i?^- = 8 — 11 or i?j = 1, 2 in dependence of the available range in the fitted direction, a sum of 
two exponentials gave always stable fits with an optimal compromise between acceptable x^/dof and statistical errors, 
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and with Vq well separated from higher excitations by a large energy gap. Subsequently, the resulting potentials^ were 
carefully fitted to eq. and the values of the best fit parameters with its errors from jackknife analyses of the data 
subsamples are listed in table IV. 



index 






A, 


Cij 


1^ _ ibir A 
9r = ^sj 




tl, Wst 
tl, Wts 
tl, Wss 


0.0335(12) 
0.0346(3) 
0.0358(7) 


0.626(67) 
0.1479(55) 
0.639(27) 


0.044(10) 
0.0401(15) 
0.0238(81) 


0.0832(4) 
0.02763(6) 
0.1105(2) 


0.561(19) 
0.5800(43) 
0.600(12) 


0.575(38) 
0.605(17) 
0.592(20) 


11, Wst 
11, Wts 
11, Wss 


0.0354(8) 
0.0351(3) 
0.0360(7) 


0.593(37) 
0.1651(41) 
0.623(29) 


0.0292(68) 
0.0372(7) 
0.0269(56) 


0.0873(4) 
0.02768(5) 
0.1111(2) 


0.592(14) 
0.5881(50) 
0.602(12) 


0.582(28) 
0.603(20) 
0.597(21) 


t2, Wst 
t2, Wts 
t2, Wss 


0.0336(2) 
0.0343(1) 
0.0345(3) 


0.594(26) 
0.1401(19) 
0.594(13) 


0.0332(65) 
0.0390(9) 
0.0346(29) 


0.0833(1) 
0.02776(2) 
0.1110(1) 


0.5622(35) 
0.5739(14) 
0.5781(54) 


0.562(15) 
0.5932(78) 
0.5781(72) 


12, Wst 
12, Wts 
12, Wss 


0.0347(2) 
0.0338(1) 
0.0344(4) 


0.555(19) 
0.1429(12) 
0.557(16) 


0.0284(56) 
0.0342(9) 
0.0303(28) 


0.0878(1) 
0.02792(2) 
0.1117(2) 


0.5821(29) 
0.5657(11) 
0.5761(69) 


0.570(14) 
0.5621(64) 
0.574(12) 


p2, Wst 
p2, Wts 
p2, Wss 


0.0339(1) 
0.0343(1) 
0.0345(3) 


0.576(13) 
0.1428(11) 
0.5810(98) 


0.0322(36) 
0.0362(2) 
0.0307(19) 


0.0851(1) 
0.02780(1) 
0.1112(1) 


0.5679(21) 
0.5742(12) 
0.5780(43) 


0.5645(92) 
0.5845(50) 
0.5756(77) 



TABLE IV. All Yukawa fit parameters of tie static potentials, calculated from space-time (ij = st and ij = ts) and 
space-space (ij = ss) Wilson loops. The renormalized coupling gjt{l / my^ij) is explained in the text. 



We only used uncorrelated fits in the present context, because the size of the Wilson loop extensions does not admit 
much variation in the fit intervals. In some cases the smallest distances i?j = 1 or i?,; = 1, 2 were omitted to have a 
satisfactory x^/dof This supports the experiences from earlier work |lj] that the lattice correction Gij may be not 
adequate enough for our data. A more thorough inspection of the fit results hints at a renormalization of = 0.5 
on the C'(15%)-level, and from the validity of Aij ~ Dij one can judge, how good the assumption of a one gauge 
boson exchange really is. The space-like potentials from Wst and Wss lead to consistent numbers, but the discrepancy 
between the screening masses mv,ij € {mv,st,'mv,ss,'mv,ts} and the gauge masses of the preceding subsection is often 
larger than expected. When comparing the two lattices, we observe only small finite-volume effects in g^, but the 
mv,ij still differ outside their — even larger — standard deviations. However, as we will see below, these effects seem 
to cancel to a great extent in the mass ratios we are mainly interested in. 

For the sake of completeness we also discuss a local definition of the renormalized gauge coupling, which goes back 
to refs. |3l|,^ and has been applied to the isotropic SU(2)-Higgs model in |l^Jl^. Since the short-distance potentials 
turn out to deviate from a pure Yukawa ansatz, we set 



IGtt aiVij{Ri) - ajVijiRi - d) 
3 lij {mv,ij ,Ri — d) — lij {mv,ij , Ri) 



(94) 



at distance Ri with mv,ij as screening masses from the large-distance fits to (pi|). Ri is the solution of the equation 



1 

Ri 



1 



lijlmy^ij, Ri — d) — Iij{mv,ij, Ri) 



(95) 



and is interpolated to the physical scale i?o,ij = l/^v.ij, giving the typical interaction range of the potential. Eq. ( p5| ) 
is motivated by requiring the force aiVij{Ri) in the continuum limit ( |88| ) to be equal to the finite difference 
[aiVij{Ri) — aiVij{Ri — d)] /d as would follow from (p9|). This improves the naive choice Ri — d/2 to tree- level p^, 
because it compensates for lattice artefacts of order 0(a^/rf). The results for d = 1 are collected in the last column 
of table IV and agree with g\ from the global definition. The errors contain the statistical errors of the potentials, the 
(ever dominating) uncertainties in the masses, and systematic errors by accounting for the sensitivity to a quadratic 
i?o,ij~interpolation with three neighbouring points instead of a linear one with only two points. 



^ More precisely, the potentials have to be rendered dimensionless before, i.e. in view of eqs. ( |86[ ) and (^) one has to attach 
a factor ai/aj. 
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Rotational symmetry now implies that the renormalized gauge coupling and mv.ij l^i should be independent of i 
and j. For g|j. this is obviously true, and in analogy to (|8^) a further kind of lattice spacing anisotropy from the ratios 
of screening masses is 

mv,st ^ mv,ss ,„„^ 
= — or = — ■ (96) 

mv.ts mv,ts 

Its values in all simulation points are quoted in table In contrast to the masses themselves, they show rather good 
consistency and are hardly affected by the finite volume. 



quantity 


tl 


11 


t2 


12 


p2 


= mv,st/mv,ts 


4.23(47) 


3.56(26) 


4.24(18) 


3.88(14) 


4.033(96) 


£,v — rnv,3s/mv,ts 


4.32(24) 


3.76(20) 


4.24(12) 


3.89(13) 


4.068(80) 


via matching 


4.250(77) 


3.923(62) 


4.179(38) 


3.915(52) 


4.028(31) 



TABLE V. Errors for the lattice spacing anisotropy obtained as ratios of the corresponding screening masses are computed 
from their jackiinife samples. The direct matching of the potentials is described in the text. 



The errors of ^h, and are relatively large. This is caused by the fact that they are determined as ratios of 
masses with individual statistical errors. The jackknife errors quoted are obtained from the jackknife samples for the 
mass ratios themselves. Calculating the errors from the mass errors using error propagation would result in even larger 
error estimates. Inspired by a method found in ref. one can obtain even smaller errors instead, if f is directly 
determined by a matching of the space- and time-like secondary quantities, without any reference to the correlation 
lengths extracted from them afterwards. We have realized this proposal for the static potentials in space (Vst) and 
time (Vts) direction. To begin with, we calculated the corresponding continuum potentials 



Vcont,ij{Ri) = Vij{Ri) ~ Cij — DijGij {mv.ij , Ri) 



(97) 



since the lattice sum lij in 



) is only meaningful for integer Ri. Constant and lattice correction terms in lowest 
order are found from eqs. ( |89| ) and ( |9ll) to be 



Cij -t- DijGij {mv.ij , Ri 
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h] {mvi] , 0) - Iij {mv.ij , Ri) + 



-rav.ijRi 



(98) 



while solely in the subtraction step gj^ and mv^ij were taken from table IV . Hence the matching condition reads 



Vcc 



t{Rs)^C-Vcont.ts{Rt/i), ^V=i- 



(99) 



It was fulfilled by fitting the space-like continuum potential to a Yukawa shape —Ae~™'^/x -I- C in imitation of (p8|), 
equating the fit function at arguments Rt/S, with the time-like potential data times a constant, and solving every 
possible equation pair for ^ and c. The final <^v-values given in the last row of table are averages over all such 
solutions along that i?i~interval, in which the two potentials have their characteristic slopes, and interchanging the roles 
of V^ont.st and Vcont,ts in eq. (|9^) always enabled a useful cross-check. As exemplarily reflected in the perturbative 
simulation parameters on the larger lattice in figure ^, the deviation between the curves then becomes uniformly 
minimal in their whole range. 

The lattice spacing anisotropics from this potential matching resemble the screening mass ratios, but the errors 
from a repetition of this procedure with 1000 normally distributed random data are indeed smaller. Moreover, is 
fully compatible with S^h and S^w in the previous subsection at the perturbative values of the coupling anisotropics. 
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FIG. 6. Matching of the (subtracted) lattice potentials at the perturbative parameters. The error bars are smaller than the 
symbols. 



D. Evaluation of the non-perturbative asymmetries and comparison with the perturbative result 

We have determined the lattice spacing anisotropics <^,; from Higgs (i — H) and gauge [i — W) boson correlation 
functions and static potentials (i = V) at different pairs of coupling anisotropy parameters. Since 7^ has been held 
fixed, each is looked upon as a function of 7/5, and the requirement of space-time symmetry restoration suggests 

the existence of a unique coupling anisotropy 7^"'''', where all possess the same value ^^"p'. This defines the 
non-perturbative anisotropy parameters. 

Therefore, we linearly interpolate the numbers = ^ii'jp.j) at the three values jpj of eq. ( |8^ ) within their errors 

A^ij to a matching point (7^"'''', C*'"^'') by minimizing the sum of squares 



EE 



(np) 



AC, 



(100) 



with respect to Ci and the common fit parameters 7^"'''' and ^("p^. We obtain the final results 

82 X 12 X 48 : 7^""^ = 3.911(43) , C'"p) = 4.040(35) 

8^ X 16 X 64 : 7^"^^ = 3.920(19) , ^("p) 4.040(26) 

with errors coming from 5000 normally distributed random data. Figure |^ illustrates that both points agree with the 
simulated Ci at the perturbative 7/3-value as well as with the perturbative point itself, and finite-size effects appear to 
be remarkably small. 



(101) 
(102) 
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lattice: 8x8x12x48, 7,= 4.0, S = 8.0, A = 0. 000178, k- = 0. 10662 

4,5 1 \ r h \ \ 

4.3 - , - 




lat tice: 8x8x16x64, y,=4.0, ^ = 8.0, A = 0.000178, /c = 0. 106 62 




3.0 I \ \ \ \ \ I 3.0 I \ \ \ \ \ I 

3.75 3.80 3.85 3.90 3.95 1.00 4.05 3.75 3.80 3.85 3.90 3.95 4.00 4.05 

coupling anisoiropy coupling anisotropy 

FIG. 7. Final ^-evaJuation for both lattices from the three simulation points, whose equal abscissas are slightly displaced 
for better visualization. The inserts show the average matching points and its error ellipses, which enclose both the numerical 
estimates at 7/3 = 3.919 and the perturbative result. 



It reimains to be mentioned that (102) includes the ^y-values — which incidentally were not available at 7^ — 3.919 
for the smaller lattice — from the matching of the potentials. Using the weighted averages of the two screening 
mass ratios in table ^ in place of the former, we get the similar results 7^"'''' = 3.921(38), ^("p' = 4.038(29), and 
7^''P) = 3.921(19), ^("P' = 4.038(26), respectively. 

All estimates signal a perfect confirmation of the perturbative results 7^'''' = 3.919 and ^'^p^ = 4.052 calculated 
in section II. and quoted in item e. at the end of subsection II. B. There is no evidence that the unknown higher- 
order corrections in and A could lead to any visible modifications, which would make the applicability of one- loop 
perturbation theory to the anisotropy parameters doubtful. In conclusion, the non-perturbative contributions can not 
be resolved within the intrinsic errors of numerical simulations, and as a consequence, the perturbative choice of the 
anisotropy parameters in investigations with the SU(2)-IIiggs model with asymmetric lattice parameters is justified 
also for Higgs masses mn > 80 GeV. 



IV. SUMMARY 



In summary, we have worked out the complete one-loop perturbation theory of the SU(2)-IIiggs model on lattices 
with asymmetric lattice spacings in gauges. We have determined the critical hopping parameter and the coupling 
asymmetries in one-loop perturbation theory, as a function of the asymmetry parameter ^. We have proven by explicit 
calculations the gauge independence of these results in gauges. We have perturbatively studied the approach to 
the continuum limit of the finite temperature theory and have determined the optimal choice of ^ ensuring the most 
economical lattice simulation for a given precision determination of the physical parameters. 

To test the relevance of the perturbative results to non-perturbative studies we have determined the non-perturbative 
coupling anisotropics using lattice simulations. Three channels have been studied, namely Higgs and W-masses as 
well as the static potential. For our parameters, i.e. Higgs mass near 80 GeV, w 0.5 and ^ = 4 the perturbative 
results agree with the non-perturbative determination within the (high) accuracy of the latter. This result opens 
the possibility to perform lattice simulation using the perturbative coupling anisotropics without the need of a non- 
perturbative determination. In particular our results are essential to study the electroweak phase transition for Higgs 
boson masses around or above 80 GeV and determine the properties of the hot electroweak plasma. 

This work was partially supported by Hungarian Science Foundation grant under Contract No. OTKA-T016248 
and OTKA-T022929 and by Hungarian Ministry of Education grant No. FKFP-0128/1997. 
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